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We present results for Drell-Yan production from the Geneva Monte-Carlo framework. We 
combine the fully-differential NNLO calculation with higher-order resummation in the 0-jettiness 
resolution variable. The resulting parton-level events are further combined with parton showering 
and hadronization provided by Pythia 8. The 0-jettiness resummation is carried out to NNLL 7 , 
which consistently incorporates all singular virtual and real NNLO corrections. It thus provides 
a natural perturbative connection between the NNLO calculation and the parton shower regime, 
including a systematic assessment of perturbative uncertainties. In this way, inclusive observables 
are correct to NNLO, up to small power corrections in the resolution cutoff. Furthermore, the 
perturbative accuracy of 0-jet-like resummation variables is significantly improved beyond the parton 
shower approximation. We provide comparisons with LHC measurements of Drell-Yan production 
at 7TeV from ATLAS, CMS, and LHCb. As already observed in e + e _ collisions, for resummation- 
sensitive observables, the agreement with data is noticeably improved by using a lower value of 
a s (M z ) = 0.1135. 


1. INTRODUCTION 

There are several different ways to obtain theoretical 
predictions for collider processes involving strongly inter¬ 
acting particles, namely fixed-order (FO) perturbation 
theory, resummed perturbation theory, and predictions 
using parton shower algorithms. 

In FO perturbation theory, the perturbative expansion 
is carried out to a given order in the strong coupling 
constant a s . The leading order (LO) corresponds to the 
lowest order in a s required to obtain the desired process. 
The next-to-leading order (NLO) then includes all terms 
of 0(a s ) relative to the LO result, and so on. 

When the process or observables in consideration in¬ 
volve large ratios of physical scales, FO perturbation the¬ 
ory can become unreliable or even break down, since for 
each order in a s one can encounter up to two powers of 
logarithms of large scale ratios. In this case, the most 
precise perturbative predictions are obtained by resum¬ 
ming the logarithmic terms to all orders in a s . At leading 
logarithmic (LL) order one includes the most dominant 
terms of order a” In 2 ™. The next-to-leading logarithmic 
(NLL) order also includes the next-largest subleading log¬ 
arithmic terms, and so on. Like at FO, the resummed 
perturbative expansion happens in a systematic fashion, 
which is described in more detail later on. 

Parton shower (PS) algorithms formally work in the 
strongly-ordered limit, where each subsequent emission 
happens at a smaller resolution scale than the previous 
one. They start from a given LO prediction and effec¬ 
tively multiply it by splitting probabilities for each addi¬ 
tional emission. This allows in principle any observable 
to be predicted at LL (and including some NLL effects) 
using a probabilistic Markov-Chain algorithm. 

Each type of prediction has its advantages as it pro¬ 


vides the most accurate description in the limit where 
other predictions are not applicable. FO predictions are 
necessary for a precise description of additional hard (i.e. 
energetic wide-angle) emissions. These are not correctly 
described even at LO by either resummation or parton 
showers, since both of these approaches rely on a kine¬ 
matic expansion in a small resolution variable. For ob¬ 
servables sensitive to many soft and collinear emissions 
generating large logarithms, FO predictions are not suit¬ 
able and resummation and parton showers are neces¬ 
sary. Resummed perturbation theory allows one to sys¬ 
tematically carry out the resummation beyond the LL 
or strongly-ordered limit and is therefore more accurate 
than parton showers. On the other hand, the higher- 
order resummation applies at the level of observables that 
are still sufficiently inclusive, i.e., that marginalize over 
many emissions and are only sensitive to a small number 
of physical scales. To obtain a prediction of the final state 
that is fully exclusive in all emissions, parton showers are 
required. In particular, parton shower predictions allow 
one to further attach a hadronization model to generate 
fully exclusive hadron-level events. These are an essen¬ 
tial requirement for experiments to be able to simulate 
their data and study detector effects, and allow for the 
most direct comparison with experimental data. 

Clearly, to obtain the best possible predictions, it is de¬ 
sirable to combine the different theoretical descriptions in 
such a way that one benefits from the advantages of each. 
In particular, this allows for reliable predictions also in 
the transition regions in between different parametric 
limits, which are often important in practice. The (often 
difficult to answer) question which parametric regime is 
the most appropriate in a given phase-space region then 
becomes much less relevant for practical purposes. 

Combining resummed and fixed-order perturbation 
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theory is a standard procedure in higher-order resummed 
calculations, where it is well understood how to systemat¬ 
ically match the pure resummed result to the fixed-order 
result in the appropriate fixed-order limit. 

The goal of combining higher-order perturbative calcu¬ 
lations with a parton shower is to obtain fully-exclusive 
hadronized events that also have the perturbative accu¬ 
racy of the higher-order calculations. So far, there have 
been many advances to combine FO calculations with 
parton showers. The combination of LO predictions for 
several multiplicities with parton showers was pioneered 
in Refs. mm- The combination of generic NLO calcu¬ 
lations with parton showers started with the methods of 
MC@NLO [3l|4] and Powheg BHZj, and is being contin¬ 
uously developed further, see e.g. Refs. mm- By now 
there are several methods available to combine multi¬ 
ple NLO calculations with parton showers H2H2Z]. More 
recently, combinations of specific Drell-Yan like NNLO 
calculations with parton showers have been presented 
in Refs. [55H32| . A general method to construct NNLO 
calculations matched with parton showers has been dis¬ 
cussed in Ref. ;33], and in which the method of Ref. [28] 
follows as a special case. 

In the above methods, the primary goal is to improve 
the perturbative accuracy of parton shower Monte Carlo 
programs for inclusive FO observables, while the de¬ 
scription of the resummation region is left to the parton 
shower. Formally, this combination amounts to matching 
the FO calculation to the LL resummation in the parton 
shower’s resolution variable. The perturbative accuracy 
criteria that should be satisfied at NNLO+PS have been 
discussed in detail in Ref. [351 . 

Geneva goes beyond this by also incorporating the 
higher-order resummation in a suitable resolution vari¬ 
able that separates the FO and parton shower regimes. 
The perturbative matching essentially happens now from 
FO to resummation to parton shower. By carrying out 
the resummation to the appropriate order consistent with 
the FO accuracy, it effectively mediates between the FO 
calculation and the parton shower and bridges the per¬ 
turbative gap between them. This provides several key 
benefits: 

• A systematic combination of NNLO calculations 
with parton showers, which extends to arbitrary 
processes. 

• An improved perturbative accuracy of (sufficiently 
inclusive) resummation observables, 

• an improved description of the matching/transition 
between resummation and fixed-order regimes, 

• a systematic assessment of perturbative uncertain¬ 
ties in the resummation and matching. 

An implementation of our approach for e + e^ collisions 
was presented in Ref. [24], where NLO calculations for 
e + e _ —)• 2 jets and e + e~ —)• 3 jets, and the NNLL' re¬ 
summation for thrust (2-jettiness) as resolution variable 


are combined with the parton showering and hadroniza- 
tion provided by Pythia 8 J34] [35]. For all the tested 
2-jet-like resummation-sensitive observables, the predic¬ 
tions from Geneva closely match the corresponding ex¬ 
act NNLL* 7 ) resummed calculations, and after hadroniza- 
tion are in excellent agreement with LEP data. Although 
not emphasized there, the results in Ref. [32j already in¬ 
clude all singular NNLO corrections, and are formally 
accurate to NNLO up to power corrections in the resolu¬ 
tion cutoff. 

In this paper, we present results for Drell-Yan produc¬ 
tion pp 4 7/Z -) £ + £~ . Our implementation combines 
the general construction of Ref. [33] with the inclusion of 
higher-order resummation in the jet resolution variable 
similar to Ref. [21]. We combine the fully-differential 
NNLO calculation for Drell-Yan with the NNLL' resum¬ 
mation for O-jettiness To (aka beam thrust 36]) with 
Pythia 8. The FO accuracy of the generated events 
is NNLO for pp —> £ + i~ , NLO for pp —> i + t~ j, and LO 
for pp —1 l + l~jj (which we refer to as NNLOo, NLOi, 
and LO 2 , respectively) while the 7o distribution is ex¬ 
actly reproduced at NNLL'. We show that by includ¬ 
ing the higher-order resummation for To as the resolu¬ 
tion variable, the perturbative accuracy for other 0-jet- 
like observables in the resummation region is significantly 
improved, and our predictions agree well with dedicated 
higher-order resummed calculations. 

Comparing our results to Drell-Yan measurements 
from ATLAS, CMS, and LHCb, obtained during the 
7TeV Run at the LHC, we obtain overall good agree¬ 
ment with the data for rapidity distributions, inclusive 
and exclusive jet cross-sections, as well as the transverse 
momentum distribution of the hardest jet. For the trans¬ 
verse momentum of the vector boson and the related 
(jf distribution between the leptons, our default settings 
with a s (Mz) = 0.118 predict a somewhat harder spec¬ 
trum. We observe a much better agreement with the data 
when using a lower value of the strong coupling constant 
a s {Mz ) = 0.1135, as was obtained from fits of N 3 LL' 
resummed calculations of event shapes to LEP data [33- 
00]. The same was also observed for resummation ob¬ 
servables in our e + e _ results in Ref. [35] . 

This paper is organized as follows. In Sec. [3] we discuss 
the theoretical framework of Geneva. In Sec. 12 Al we 
provide an overview of the general construction and the 
master formulas for the partonic jet cross sections that 
serve as the main perturbative input, and in Sec. |2 B| 
we discuss specific choices and implementation aspects 
that are of relevance in more detail. In Sec. [3] we de¬ 
scribe how the partonic events with up to 2 partons are 
interfaced with a parton shower in a way that avoids dou¬ 
ble counting and such that the parton shower does not 
destroy the perturbative accuracy of the partonic calcu¬ 
lation. In Sec. [i] we compare Geneva with dedicated 
higher-order perturbative calculations, and in Sec. [5] we 
compare Geneva’s results with measurements from the 
LHC. We conclude in Sec. [G] 
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2. THEORETICAL FRAMEWORK hand notation 


A. General Setup 


d$M 

dd>jv 


d$M - $jv(®m)] • 


(3) 


In this subsection we discuss the general construction 
we use. Our discussion here follows that of Ref. [32) and 
we refer to that paper for more details. 

Drell-Yan production at NNLO receives contributions 
from partonic processes with up to two final-state par- 
tons. In a standard FO calculation, the phase space for 
these partonic contributions is integrated over in such a 
way that the cancellation of virtual and real IR diver¬ 
gences happens at the level of observables. 

In contrast, an event generator is meant to produce 
physical events, which means that all IR divergences 
should cancel on a per-event basis. This implies that an 
IV-parton event should correspond to an IR-finite par¬ 
tonic IV-jet cross section, which is fully-differential in the 
corresponding partonic IV-jet phase space. That is, each 
generated event represents a point in an IV-jet phase 
space, rather than an IV-parton phase space, and each 
four-momentum in the event provides the energy and di¬ 
rection of a partonic jet. 

Therefore, the basis of the Geneva MC framework 
is the formulation of the perturbative inputs in terms 
of “Monte Carlo (MC) cross sections”, which are well- 
defined partonic jet cross sections according to which the 
events are distributed, and which can be systematically 
computed to the desired perturbative accuracy in FO and 
resummed perturbation theory HH HU Hg. 

The MC cross sections are defined in terms of an IV-jet 
resolution variable Tn and formally include the contri¬ 
butions of an arbitrary number of unresolved emissions 
below a resolution cutoff Tn < T£ ut - In the present case, 
we require events with 0, 1, and 2 partons, which are 
distributed according to 

Hrr MC 

$ 0 events: ^(7^ cut ), 

drr MC 

events: (To > T 0 cut ; 7) c:ut ), 

der MC 

$ 2 events: —p (T 0 > T 0 cut , T > T^) ■ (1) 

<W2 


The exclusive 0-jet MC cross section is defined by 7o < 
7() CUt , the exclusive 1-jet MC cross section by 7o > 7)J ut 
and 71 < T{ ut , and the inclusive 2-jet MC cross section 
by 7o > Ty ut and T\ > 7] cut . In this way all of partonic 
phase space is covered. Adding the 1-jet and 2-jet events, 
we can also define the inclusive 1-jet MC cross section as 


d<T> c j 

d$i 


(To > T 0 cut ) 


daf c 

d$i 


(To > T 0 cut -, 


T cut ) 


( 2 ) 



d$2 dugi 

d$i d$ 2 


(r 0 >7? ut ,r 1 >rr t ), 


which is defined by 7o > 7f ut and does not depend any¬ 
more on Tj cut . In Eq. ([2]) we have made use of the short- 


Since the partonic jets are represented by on-shell 
quarks and gluons in the partonic events, the par¬ 
tonic jet definition used here is quite different from an 
ordinary jet algorithm. It depends on a phase space map 
‘FaK^’m) that projects the M -body phase space of unre¬ 
solved emissions onto the distributed "Fat points, where 
N < M. This phase space map must be IR safe, such 
that the MC cross sections are IR finite. For example, at 
fixed NNLO the exclusive 0-jet cross section is given by 

drr MC 

^-(7^ cut ) = (So + K, + ^o)($o) 

+ J^ Q (Bi + V 1 m)e[T 0 ^ 1 )<T 0 c ' lt } 

+ j ^ s 2 ($ 2 ) 0[r o ($ 2 ) < r 0 cut ], (4) 

where Bn contains the IV-parton tree-level contributions, 
Vn the IV-parton one-loop contributions, and Wo the 
two-loop contribution. There is considerable freedom in 
the precise definitions of Tn and the &n($m) map, and 
we will choose particular definitions based on their theo¬ 
retical properties and simplicity. 

Using the events in Eq. 0 . the cross section for any 
observable X is given by 

r Hrr M0 

a(X) = y d $o^(7^ cut )M x ($„) (5) 

r Hrr MC 

+ J d$! ^(To > 7J ut ; Tj cut ) M x (*i) 

I" dcr MC 

+ J d$2 > T 0 CU \ T > Tn M x (T 2 ) , 

where Mx(^n) is the measurement function that com¬ 
putes the observable X for the IV-parton final state <Fat. 
This cross section is not identical to the exact fixed-order 
result, because for any unresolved emissions the observ¬ 
able is calculated on the projected phase space point 
rather than the exact <F M . This is a funda¬ 
mental limitation inherent to generating IR-finite events 
beyond LO and for a more detailed discussion we refer to 
Ref. [33] • The difference vanishes in the limit 7)^ ut —> 0. 
Hence, we want to choose 7]^ ut as small as possible to 
have a maximally exclusive description. For small 7]^ ut , 
however, the MC cross sections contain large logarithms 
of Tn and 7 ~N Ut i requiring resummation to obtain phys¬ 
ically meaningful predictions for them. This is precisely 
what a parton shower would do at LL (also including 
some NLL effects), in which case T^ ut would play the 
role of the shower cutoff. In Geneva we are improving 
over this by resumming these logarithms at higher accu¬ 
racy. 
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1. 0/1-jet separation 


To discuss the construction of the MC cross sections, 
we first consider the separation between 0 and 1 jets 
which is determined by the 0 -jet resolution variable 7o- 
We write the exclusive 0-jet and inclusive 1 -jet MC cross 
sections as [24l l33l 


dg“° 

d $ 0 


(T 0 cut ) 


da^ 

d$i 


(To > T 0 cut ) 


da 5 esum 

d$ n 


(V 


da " 118 


match 


d$n 


j nons 

uu 0 f'-rcut \ 

d$ n (/ ° 


da| e f um 

d$i 

da |7 


9(T > r 0 cut ) 


match 


d$i 


-(To >r 0 cut ) 


XU >1 

d$i 


(To > T 0 cut ). 


-(7? ut ) 


( 6 ) 


Here, daQ esum is the spectrum differential in $o with 
the dependence on 7jj cut resummed to a given logarith¬ 
mic accuracy, while the remaining two terms give the 
matching corrections required to reproduce the desired 
FO accuracy. Their precise form depends on the fixed- 
order content of the resummed contribution. The singu¬ 
lar matching dag lng match contains all contributions that 
do not vanish as T^ ut —> 0 , while the nonsingular match¬ 
ing daQ° ns only contains contributions that vanish in this 
limit. The analogous separation is done for da|°, where 
da>f um resums the differential dependence on %, the 
singular matching contains contributions that diverge at 
least like l/7o for To —s► 0 , and the nonsingular matching 
contains contributions that contain at most integrable 
singularities for To —>■ 0 . 

We carry out the resummation of the To dependence 
to NNLL'. Thus, the resummed contributions are given 
by 


da 5 esum 


d$o 


(T 0 cut ) 


d(T resu m 


d$i 


da NNLL ' 

d $ 0 

da NNLL ' 

d$ 0 d7o 


To' 


T($i). 


(7) 


Here, da NNLL /d7od<l>o is the resummed differential To 
spectrum and da NNLL /d<i>o(7J ut ) its cumulative inte¬ 
gral, see Sec. |2 B 2| below. Note that the two-loop vir¬ 
tual corrections, which at FO are proportional to 5(To ), 
are fully incorporated at NNLL'. This means they are 
properly spread to nonzero values of To and contribute 
to the differential spectrum as dictated by the resum¬ 
mation. The To resummation is naturally differential in 
$o, so we can directly use the resummed cumulant in 
dp-resiim. qy, ma k e the 7 ^ spectrum fully differential in 
5q, we have defined a normalized splitting probability 
T(<Iq) which satisfies 


" d$i 

d<J>odT) 


T($i) = l. 


( 8 ) 


Since To can be considered as part of the radiation phase 
space d<f>i/d<i>o, the integration here is effectively over 
two remaining radiation variables, e.g. an energy split¬ 
ting and an azimuthal angle. Thus, the integral of the 
splitting probability over all <Eq points restricted to any 
fixed values of $o(3q) and TX'&i) is equal to 1. This will 
be discussed in more detail in Sec. I2B4I 

At NNLL', the resummation fully incorporates all sin¬ 
gular contributions to O(al), implying that the singular 
matching vanishes, 


d^>r 


i sing match 

_ fq-CUt\ _ Q 

d$o (/ ° } ’ 

match 


d$i 


-(To >7? ut ) = 0. 


(9) 


At our desired NNLO accuracy, da|° and da|° must 
be correct to NNLOo and NLOi, respectively, which de¬ 
termines the nonsingular matching corrections to be 


da£ ons 

d$ n 


(To 


da. 


NNLOo 
o 


d$ D 


-(^ cut ) - 


cla^j 118 


d$i 


(To > T 0 cut ) 


da^ Gl 

d$i 


da£ esum 

d $ 0 


(T 0 cut ) 


J NNLOo 


( 10 ) 


(To > T 0 cut ) - 


da^sum 


d$i 


1 J NLOi 


9(% > T 0 cut ). 


The terms in square brackets are the FO expansions to oil 
of the resummed cumulant and spectrum in Eq. (7). The 
NNLOo result for the cumulant is given in Eq. (4). The 
NLOi result for the fully-differential spectrum is given 

by 


^ NLOi 

^ r (T 0 > T 0 cut ) = (Bi + HrXdq) 9[%Ti) > 7^1 

+ j^B 2 (<S> 2 )e[To(<S>2)>To c T- 

( 11 ) 


It depends on the projection 

^ = d$ 2 «5[$i-$r($ 2 )]e r ($ 2 ), (12) 

where 0 T ($ 2 ) defines the region of $2 that can be pro¬ 
jected onto massless dq via the IR-safe phase space map 
$i~(^ 2 )- Only this projectable region of <J > 2 is included 
in da>° ns /d3q, while the remainder will be included in 
the $2 events below. 

The fact that the NNLL' resummation reproduces all 
singular terms of the full FO result to O(al) implies 
that the resummed expanded result in da>° 1 ns /d4>i acts 
precisely as a differential NNLO 7o-subtraction 113] (see 
also Ref. M)- The cancellation of the singular terms 
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between the FO and resummed expanded results is not 
point-by-point in $1 but nonlocal, which means it only 
holds “on average” upon integrating over d<l>i/d$od7o- 
As discussed in detail in Ref. [15] , the cancellation can be 
made more local by considering more-differential resum¬ 
mations (e.g. utilizing the results of Refs. [45ll47j ). For 
the cancellation to be point-by-point in 7o, d< 7 >^° 1 /d < f’i 
has to reproduce the right singular 7o-dependence when 
projected onto d7od<ho, which means the map ($ 2 ) 
has to preserve the value of 7o, 

T 0 ($T(<f> 2 )) =To($ 2 ). (13) 

In other words, if this was not the case, we would not ob¬ 
tain the correct nonsingular corrections for the 7o spec¬ 
trum, and this would destroy its NNLL' accuracy. The 


underlying reason is that we perform the FO calcula¬ 
tion in an event generator form as in Eq. ([ 5 ]), and so to 
compute the FO To spectrum we always project some 
$2 points, namely those below 71 < 71 cut , onto $1 and 
only then onto To- More details on this map are given 
in Sec. 12 B 61 


2. 1/2-jet separation 

Next, we consider the separation of the inclusive 1- 
jet cross section into an exclusive 1-jet MC cross section 
and an inclusive 2-jet MC cross section, using the 1-jet 
resolution variable 71. We write them as 


dcr f 0 r7 - ^ -T-CUt. -T"CUt\ dcr i 

— (To>T 0 ,7j ) =- 


da%% 

dd>2 


(To > 71 cut ,7! > 71 cut ) = 


d$i 

do">| um 


-(To>T 0 cut ;7r t ) + 


^£j.match 

d<E>i 


(To>T 0 cut ;Tr t ), 


(To > r 0 cut ) 0(7i > rr‘) + 


d(7™f tch 


(To > T 0 CU \71 > 71 cut ) 


do-resum and d(T ™sum 


d$ 2 U 7 ^ 1 ' d$ 2 

contain the resummation of the 7j cut and differential 71 dependence. The der 


( 14 ) 


match and 


Here 

d^match corda j n the matching corrections required to achieve the desired FO accuracy. 

For the 71 resummation we limit ourselves to LL, which is sufficient for matching to the parton shower. At this 
level, we can write the resummed contributions as 


do-f 1 


d$i 


-(7o > 71 cut ;71 cut ) = 


j—resum 
U(J >2 

d <& 2 


(To > T 0 mt ) = 


dcr >i 

d$i 

dofi 

d<f>i 


Ui($i, 71 cut )0(T o >71 cut ), 


U[^ 1 ,T 1 )9(T 0 >T 0 cut ) 


$i=$7 ($ 2 ) 


7>($ 2 )0(7! >71 cut ). 


( 15 ) 


Here, [/i($i, 71 cut ) denotes the Sudakov factor that resums the dependence on 71 cut to at least LL accuracy, and 
f7{($i,7!) denotes its derivative with respect to 71 cut and resums the differential 71 dependence. Its details are given 
in Sec. 12 B 51 

The differential 71 resummation in der> 2 is evaluated at the projected $1 point (H> 2 )- As we will see below, this 
is required so the To resummation is evaluated at the correct To value. The normalized splitting probability V( $> 2 ) is 


2 B 4 


defined analogous to 7 ? ($ 1 ) in Eq. ([8]) but for —>• 4> 2 and using the $£($ 2 ) m &p, with the details given in Sec.^ _ 

The 71 resummation acts on top of the overall da?, 1 /d^i in Eq. p~5| , which is the inclusive 1-jet cross section in 
the 71-singular limit. Its NLOi expansion therefore has to be 

dc>i 


d4>i 


= (B 1 


J NLOi 


Hr)($r) + J^ C 2 ($ 2 ) = (B 1 + F L c ')(4 , i). 


( 16 ) 


(Its full form will include the To resummation and is derived below.) Here, C 2 (& 2 ) denotes a standard NLO sub¬ 
traction HHi,:®] that reproduces the pointwise singular behavior of B 2 ($ 2 ). We use the FKS subtractions [H] in our 
implementation. Note that the subtraction term here uses its own projection d$ 2 /d < t>f = d<E> 2 ($ 2 )]. Using 

Eqs. (15) and (16) and requiring that dnf 0 and der >2 are correct to NLOi and L0 2 , respectively, then determines the 
matching corrections to be 


^0.match 


j —match 
aa >2 

d$2 


itch [' Hd) rid) 

d(To > To cut ; 71 cut ) = J H 2 (4>2 ) 0[To(4> 2 ) > T 0 mt ] 0(T[ < 71 cut ) - ^ C 2 ($ 2 ) 9(% > T 0 cut ) 

- B^i) t/i (1) ($!, 71 cut ) 9(T 0 > T 0 cut ), 

) = {H 2 ($ 2 ) [1 - © r ($ 2 ) 0(71 l < 71 cut )] 

-B 1 ^l)U { 1 1) \<^'[,T 1 )V^ 2 )9(T 1 >Tr t )}0[To^ 2 )>T o cut }. 


( 17 ) 


(To > T 0 CU \T > T/ 


( 18 ) 
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Here, (4>i, 71 cut ) denotes the(7(a s ) term in the expansion of the t/ 1 ^($i,Tj cut ) Sudakov factor (or its derivative). 

These terms cancel the leading double-logarithmic terms in 71 in the FO pieces (the B 2 contributions). As long as we 
only include the LL 71 resummation, there will be some subleading single-logarithmic terms present in the matching 
corrections. 

So far, we have discussed the construction of an additive NLO 1 +LL 1 matching. The important point is now that 
this should also include the 7o resummation which is important at small 7o- The condition for this is that the integral 
of NLO 1 +LL 1 result with the $^($ 2 ) map must reproduce the 7o-resummed result for the inclusive 1-jet MC cross 
section in Eq. That is, 


d<r>j 

d$i 


(To > T 0 cut ) = -=MT 0 > 7? ut ;7? 


d$_ 

do fi 

d$i 


■> + / 


d$2 


do">2 


(To > r 0 ci 


d$f d $ 2 

0(T O > T 0 cut ) + J ^ 2 ) 0[To($ 2 


fc ,Ti >rr‘) 


> T 0 cut ] 


^C 2 ($ 2 )0(T o >V 


(19) 


where we have used the identity 


[/,($!, 7? ut ) + J ^U[(<f>i,Ti)V($ 2 )e(Ti 


> 71 cut ) = 1. 

( 20 ) 


Inserting the expression for dcr“° in terms of dcr>f um 
and dcr"° ns in Eqs. ([7]) and (10), we obtain the result 
for dcr^/d'f’i beyond NLOi, 


dqfi 

d$i 


d<j>® um 

d$i 


+ (H 1 +E 1 C )($ 1 )- 


d^resum 

d$i 


J NLOi 
( 21 ) 


It contains the full differential 7o resummation via 
dcr>® um . The difference to the full inclusive 1-jet cross 
section is that the 7 o-nonsingular terms are now evalu¬ 
ated in the 71-singular approximation, and this mismatch 
is compensated by the 71 matching contributions. 

The above expressions completely define the fully- 
differential jet cross sections. In the next subsection we 
provide some additional details on the specific implemen¬ 
tation of the various pieces in Geneva. 


B. Implementation Details 

1. Choice of the jet resolution variables 

We choose IV-jettiness [50] as our JV-jet resolution vari¬ 
ables. It is defined as 

Tn = ^minjg a ‘PkAb ‘PkAi • Pk,---AN , ( 22 ) 

k 

where the sum over k runs over all final-state particles, 
excluding the vector boson and all its decay products. 
We use a canonical geometric measure, where qi = rii = 
(1 ,n,;) are lightlike reference vectors along the jet and 
beam directions. 

iV-jettiness is a global event shape that is explicitly de¬ 
signed as an A^-jet resolution variable, as it measures the 


r 


degree to which the final state is IV-jet-like for a given N 
and automatically clusters the final state into N jet and 
beam regions in an IR-safe way and without any depen¬ 
dence on an additional jet clustering algorithm. Further¬ 
more, it is theoretically well understood and its all-orders 
singular structure and resummation are known. 

For Drell-Yan, we need 0-jettiness 7o and 1-jettiness 
71, and we always evaluate Eq. (22) in the longitudinally 


boosted frame where the vector boson has zero rapidity. 
The reference vectors for the beam directions are always 
along the beam directions, so n a = z and fib = —z. For 
71, we also need to define a jet direction nj. For our 
calculation we only need to define it for a <I >2 event, in 
which case we can choose n j by minimizing 71. We define 
a clustering metric 


d(p ) = \P\~\P Z \, 

d(pi,P2) = \pi\ + \P2\-\pi+i>2\, (23) 

and find the smallest of d(pi), d(p 2 ), and d(pi,p 2 ), where 
Pi and p 2 are the two parton momenta in the $2 event. If 
d(pi) is smallest, we can think of pi being clustered with 
one of the beams, and then nj = Pj/\Pj\ with j 7 ^ i. If 
d{pi,P 2 ) is the smallest, we can think of p\ and p 2 being 
clustered together, and then ftj = (pi + P 2 )/\pi + P 2 \■ It 
is not hard to show that this minimizes 71 and its final 
value is determined by the minimum metric as 


2 


71($ 2 ) = - \pl\) +min{d(pi),d(p 2 ),d(pi,p 2 )} • 

k =1 

(24) 


2. The 7o spectrum at NNLL' from SCET 


The all-orders parton-level factorization theorem for 
0 -jettiness (aka beam thrust) is given by JTTi (FT] 


da SCET 

d$ 0 d71 


^—. da'-, „ f 

>) dt a dt b Bi(t a ,x a ,p:) 

ij ° 
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x Bj(tb,x b ,n) s(t 0 ~ ta ■ ( 25 ) 

Here, dcr,^/d$o is the Born cross section for the ij —► 
Z/ 7 * —> £ + £~ hard scattering. The hard function Hij{Q) 
contains the corresponding Born and virtual squared ma¬ 
trix elements, and the sum runs over all possible qq 
pairs ij = {uu,uu,dd,dd ,..The Bi(t,x) are inclu¬ 
sive (anti)quark beam functions I■161 . with the momen¬ 
tum fractions x a ,b given in terms of the total rapidity Y 
and invariant mass Q = mg+e- of the vector-boson final 
state by 


x a = 


Q_ 

E 

-‘—'cvn 


„Y 


X b = 


E, 


Q -Y 
e 


(26) 


The beam functions are computed perturbatively in 
terms of standard PDFs fj , schematically Bi = JT X () - <g) 
fj. Finally, S(k ) is the quark hemisphere soft function 
for beam thrust. 


Equation (25) is derived using Soft Collinear Effective 
Theory (SCET) fo!Mj(T| . Its key feature is that each of H , 
B , S only depend on a single characteristic scale. There¬ 
fore, there are no large logarithms in their perturbative 
expansion when each is evaluated at its own characteris¬ 
tic scale, which is given by 


Mb = Q , Mb = \J QTo , Ms = 7o ■ 


(27) 


In Eq. (251, all ingredients have to be evaluated at a com¬ 


mon scale m, which is arbitrary and whose dependence ex¬ 
actly cancels between the different functions. The renor¬ 
malization group evolution in the effective theory is then 
used to evolve each function from its own scale to /i- This 
gives the resummed To spectrum used in Eq. |7]) as 


do - 


NNLL' 


der 


d$od7o 


= ^2 £Mmb,m) 


x [Bi{x a ^ B ) <8 t/ B (Ms,M)] 
x [Bj(x b ,HB) <8 Hb(mb,m)] 
<8 [S(/j. s ) <8 ^ s ( ms , m )] , 


(28) 


where we have now written the convolutions between the 
different functions in schematic form. 

The RGE factors U_ y(mx,m) resum the large logarith¬ 
mic terms, and the functions evaluated at their own scale 
provide the boundary conditions for the evolution. At 
NNLL' accuracy, we need the boundary conditions at 2- 
loop order, and the evolution to 3(2)-loops in the cusp 
(noncusp) anomalous dimensions. All required expres¬ 
sions are known in the literature 03 UZHIS] and we do 
not reproduce them here. 


In practice, the canonical scales in Eq. (27) are ap¬ 


propriate in the resummation region, where the singular 
corrections dominate. In Fig. [l] we compare the singu¬ 
lar, nonsingular, and full results for the To distribution 
at NNLO 0 - We can see that up to 7o < 50GeV the sin¬ 
gular dominate and the nonsingular are suppressed by an 



FIG. 1. Comparison of the full, singular, and nonsingular To 
distribution at NNLOq. 


order of magnitude or more. On the other hand, in the 
FO region for larger values of 7o ~ Q, the resummation 
must be turned off since here the singular terms being 
resummed become meaningless and there are large can¬ 
cellations between the singular and nonsingular terms, 
which must be preserved to reproduce the correct FO re¬ 
sult. This is clearly visible in Fig. [l] for To 80GeV, 
where the magnitude of the singular and nonsingular are 
larger than the full result (the singular turn negative at 
the dip at To — 70GeV). For this reason the resum¬ 
mation must be turned off in this region. This is done 
by taking all scales to be equal to the common FO scale 
Ms = Ms = Mb = Mfo- A smooth transition between 
the canonical and FO limits is achieved by using profile 
scales 133IM1, where mb and Ms are smooth functions of 
To, interpolating between the exact canonical scaling in 
Eq. (271 and the FO scale. The uncertainties from the re¬ 
summation and the transition between resummation and 
FO can be estimated by choosing different sets of pro¬ 
file scales, which provide a sensible variation around the 
central scale choice. In other words, for each event we 
perform the calculation of its MC cross section with a 
different set of profile scales, producing different weights 
for the event, and the variation in these define our event- 
by-event uncertainties as described next. 


Profile scales and their variations are used in many 
SCET calculations and have been shown to provide re¬ 
liable perturbative uncertainty estimates for resummed 
and matched predictions for observables in a variety of 
contexts, see e.g. Refs. [33 EZHZI]. For our central pro¬ 
files we take 


Mb' = Mfo , 

Ms(7o) = MFo/run(7o/Q) , 

Mb( 7o) = MFoV/run(To/Q), (29) 

where the common profile function f rU n{x) is as in 
















Ref. [55] , 


/run(^) 


'Xq [l + (&/^o) 2 /4] 

X 

I (2-x 2 -x 3 )(x-xi) 2 
' _r 2(x 2 —xi)(x 3 —xi) 
1 _ (2-xi-x 2 )(x-x 3 ) 2 
2(x 3 —xi)(x 3 —x 2 ) 

1 


x < 2xo , 

2x’o < x < X \, 

Si < x < s 2 , , 

X 2 < x < x 3 , 

X 3 < X . 

(30) 


which has strict canonical scaling below x 3 and the re¬ 
summation is turned off above x 3 . From Fig.[l]we choose 
for our central scale the parameters 


Mfo = Q, x 0 = 2.5G eV/Q, 

{xi,S2,x 3 } = {0.2,0.5,0.8} . (31) 


For the FO variations we vary ufo up and down to 2 Q 
and Q/ 2, taking the maximal absolute deviation in the 
results from the central value as the FO uncertainty. For 
our resummation uncertainties, we use the profile vari¬ 
ations designed for 7”-like resummation variables from 
Ref. [70j . which give us four independent up and down 
variations for hb and ns that probe deviations from the 
canonical scaling but never violate the parametric scal¬ 
ing n% ~ and turn off beyond x 3 . In addition, 

we include two more profiles where we vary all Xi transi¬ 
tion points by ±0.05, giving us a total of 6 resummation 
profile variations. We then take the maximal absolute 
deviation in the result from the central value among all 6 
profiles as the resummation uncertainty. The total per¬ 
turbative uncertainty is obtained by adding the FO and 
resummation uncertainties in quadrature. Note that the 
profile scales depend on the underlying 4>o point through 
the value of Q, but we choose them to be independent of 
Y. 

The same resummation formalism can be used to ob¬ 
tain an expression for the resummed cumulant of the To 
spectrum in Eq. 0 


da NNLL ' 

d$ 0 


(T 0 cut ) 



da NNLL ' 

d<F 0 d7[) 


(32) 


Here, we first integrate the resummed expression for the 
To distribution in Eq. ( [28]) , and then choose the scales 
using the same profile scales but as a function of 7[) cut . 
In particular, in the resummation region, the canonical 
values are 


m = Q, hb = \[q vs = T 0 cut . (33) 

Since the operations of integrating the factorization the¬ 
orem and choosing the scales do not commute with one 
another, the resulting expression for the cumulant is not 
exactly equal to the integral of the To spectrum, where 
the scales would be chosen before performing the inte¬ 
gration. The resulting differences are formally always of 
higher order, beyond NNLL' in our case (see e.g. Ref. m 


for a detailed discussion). However, they can be numer¬ 
ically relevant, especially if one aims at preserving the 
NNLO 0 total cross section as accurately as possible. 

In principle, one could enforce the equality between the 
cumulant and the integral of the spectrum in Eq. (32) ex¬ 


actly, by simply defining the To spectrum as the deriva¬ 
tive of the 7[) CUt cumulant. However, this can give rise to 
unreliable uncertainty estimates in the spectrum as well 
as a poor description of the region of To where the resum¬ 
mation is turned off and FO predictions are valid 1371 . In 
Geneva we circumvent these problems by adding to the 
To spectrum the quantity 


=(T 0 ) 


d da NNLL ' 
dT 0 d$ 0 


(7o,/Pi) - 


der 


NNLL' 


d-F 0 d7() 


iVh) 


. (34) 


which is explicitly beyond NNLL' order. The profile 
scales nh are chosen to turn off the resummation some¬ 
what earlier than the profile scales for the rest of the 
resummed calculation, in order to preserve the accurate 
spectrum prediction in the tail of the distribution. Fur¬ 
thermore, nh and its derivative are smooth functions of 
To, which is needed because the difference of resummed 
terms in Eq. (34) is proportional to dfih/dTo [72] . The 
function k(To) is a function that smoothly interpolates 
from a constant value for To <C Q to zero at the point 
where nh{To) = Q , and further reduces the effects of the 
higher order correction terms in the tail of the distribu¬ 
tion. By adjusting the constant value at small To within 
0(1) factors, the overall size of these terms can be ad¬ 
justed (in practice this constant value is approximately 
2 ). 

These higher-order terms may be used to capture the 
correct FO scale uncertainties in the inclusive cross sec¬ 
tion upon integrating over the full spectrum. We use 
the same profile for each of the scale variation described 
above, but scaled by the value of ufo/Q (so 1 for the 
central scale and resummation scale variations, 2 for the 
up FO scale variation, and 1/2 for the down FO scale 
variation). By adjusting the constant value of k(To —1 0) 
separately for each n-Fo/Q value, we can correctly obtain 
the total integrated FO cross section for each FO scale 
variation. 

Technically, each of the resummation scale variations 
should yield strictly the same total cross section, such 
that the resummation uncertainties vanish for inclusive 
observables. However, in practice, the profile scale vari¬ 
ations are primarily meant and designed to capture the 
resummation uncertainties in the peak and transition re¬ 
gions. Implementing a constraint that they all reproduce 
exactly the same numerical integral of the resunnned 
spectrum as the central profile is very challenging, es¬ 
pecially since the resummation uncertainties in the peak 
region can be sizeable, while the inclusive cross section 
uncertainties are very small. (Alternatively, we could 
separately finetune the k(To) function for each of the dif¬ 
ferent resummation profiles.) Instead, what we will do is 
to simply drop the resummation uncertainties whenever 
we consider fully inclusive observables. 
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FIG. 2. The nonsingular cumulant for O-jettiness, as a func¬ 
tion of 7o cut . The black (blue) curve shows the NLO (pure 
NNLO) nonsingular cumulant, coming from the /i and /2 
terms in Eq. (351. 


3. NNLO nonsingular corrections 


As discussed in Sec. |2 A 1[ the fact that we are resum¬ 
ming the 7o spectrum and cumulant to NNLL' accuracy 
implies that there are no singular terms in the matching 
to the NNLO calculation. Therefore, Geneva directly 
computes the fully-differential inclusive 0-jet cross sec¬ 
tion at NNLO, up to power corrections in 7)) cut . The 
power-suppressed missing terms are precisely the 0-jet 
nonsingular contributions in Eq. (101, which we can write 
as 


der£ ons 

d<£>o 


(T 0 cut ) = h/i(7T\ ®o) + «s/ 2 (T 0 cut , $o)]77’ 


s(T 0 cut ) = 


f dCT S° nS 


(T 0 cut ). 


(35) 


The functions /fc(7)j CUt , To) can still have divergences as 
7jj cut these divergences are at most logarithmic. 

Thus, these corrections vanish in the limit 7jj cut -A 0, and 
for small enough 7)) cut they can be neglected, allowing 
to obtain NNLO accuracy without doing the full NNLO 
calculation (which would determine the function / 2 (d>o)). 

We show the size of the total correction terms 
S n ons(7() CUt ) in Fig.[2]for the NLO and pure NNLO contri¬ 
bution. At our default value of Tp ut = 1 GeV these terms 
are very small but also not completely negligible. In the 
current implementation we include the full NLO term, 
i.e., the /i($o, 7jj cut ) piece in dcrg ons (7(> cut ) by perform¬ 
ing the corresponding 0-jet NLO calculation on-tlie-fly. 
For the NNLO term, we neglect the $0 dependence and 
only include its total integral. Technically, this is done 
by obtaining the pure NNLO contribution to E nons (7jj cut ) 
from the known dependence of the total cross section on 
To ut (i.e. analogous to how we obtain Fig. [ 2 ]). This cor¬ 
rection is then included through a simple rescaling of 

d^ c /d$ 0 (T 0 cut ). 


4- The normalized splitting probabilities 


In Eq. Q we introduced a normalized splitting prob¬ 
ability 7 , ( < i>i), which is used to make the resummed 7o 
spectrum fully-differential in $1 in Eq. ( |7|. A similar 
splitting probability V(<& 2 ) is used in Eq. (|15|) to make 
the Ti resummed term fully-differential in $ 2 - As men¬ 
tioned before, the splitting probabilities are normalized, 


d<f> 


A + l 


d$7vdT/v 


^(^iv+i) — 1 ■ 


(36) 


and provide the dependence on the two extra emission 
variables required besides Tn to define a full splitting 
from —> 'Fjv+i- The two extra variables are chosen 

as usual as an energy ratio z and an azimuthal angle 
0. In the collinear and soft limit, the z variable denotes 
the ratio z = Em/E$ for a splitting M —► DS, where 
a mother particle M splits into a daughter particle D 
and a sister particle S. For ISR splittings, the daughter 
particle is the final-state particle, while for FSR splittings 
it is always the gluon for q -A qg splittings, the quark for 
g —> qq ones and the softer gluon for g —>■ gg splittings. 

We first need to decide which two particles in $jv+i 
are interpreted as arising from the splitting and which 
are then labeled with D and S. Here we have to choose 
the two particles which set the value of Tn, which are the 
two particles that are closest according to the IV-jettiness 
clustering metric, as described in Eq. (231. 


The normalized splitting probability is then given by 


V^N+l) = 


P S p{zA) 


d$ArdT/vdzd</> 




d$/> 


(37) 


where the splitting function p sp (z, </>) depends on the cho¬ 
sen splitting sp which as explained above depends on the 
phase-space point sp = sp($jv+i). With this definition, 
the splitting probability is indeed normalized: 


dd>Ar+i 


T(^n+ 1 ) 


d<f>7vdT/v 

= [ dzd ^ 


Psp{ z A) 


sp 

= 1. 


Es P fzy z p°p (*» 


(38) 


The splitting functions for FSR and ISR are chosen to 


be 


PspeFSR(z, 4>) = AP S p(z, </>) , 

PspSlSR^, 4>) = APs p("’ ’ ( 39 ) 

with AP sp (z, </>) denoting the unregularized Altarelli- 
Parisi splitting function, and f(x, g) the standard parton 
distributions function evaluated at the momentum frac¬ 
tion x and the factorization scale p. The variables Xm 
and x$ denote the momentum fractions of the mother 
and sister particles in the ISR splitting M —> DS. 
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5. 71 resummation 


The Sudakov factor di(T[ max ,7i) in Eq. (15) resums 
the 71 dependence to LL accuracy. We use the expression 
that is obtained from the 71 factorization theorem [50 
analogous to the one given for 7o in Eq. (251, where all 
ingredients are calculated at tree level, and the running 
is performed at LL only. The 0(a s ) expansion of this 
Sudakov factor has the simple expression 


tf 1 < 1) (77 nax ,Ti) 


a s (7? nax )(2CV + C A ) 

2tt 



max 


(40) 


usual FKS maps are designed to be invertible, the invert- 
ibility turns out to be a very nontrivial requirement on 
the d*/"(d> 2 ) map. 

Formally, we only need to preserve 7o for the singular 
terms, i.e., we can use any map which preserves a vari¬ 
ant of To with the same singular structure to the order we 
want to carry out the resummation, while any nonsingu¬ 
lar difference is captured by the nonsingular corrections. 
Therefore, to enable us to invert the map, we hold a re¬ 
cursive definition of To fixed, which is effectively defined 
by the action of the map, 

r 0 FR ($ 2 ) = r 0 [$n$2)]. (4i) 


The above minimal LL resummation is sufficient for 
our current purposes, where we focus on 0-jet resum¬ 
mation variables and inclusive observables and where we 
primarily need it to have a proper matching to the parton 
shower. A more precise description of 1-jet resummation 
variables will require an improved 71 resummation, which 
we leave for future work. Extending the 7i-resunnnation 
beyond LL accuracy is of course possible. One can also 
consider including the resummation through a Powheg- 
like Sudakov factor, where the ratio of the real matrix el¬ 
ement i ?2 over the born B\ would be exponentiated. We 
have verified that our underlying NLOi+LLi matched 
calculation reasonably agrees with a Powheg implemen¬ 
tation of V + 1 jet [73] in regions where the 0-jet resum¬ 
mation is not relevant. 


6. NLOi calculation and phase-space map 

The last point addressed in this section concerns the 
phase-space map <F/($ 2 ) that is used to project a par- 
tonic phase-space point $2 onto the $1 point in the NLOi 
calculation, see Eqs. EL and It is also used in the 
71 resummation in Eq. HI- Just requiring to have some 
IR-safe map as usual for fixed-order calculations (like for 
example the FKS map ($ 2 ) used in the C 2 subtraction 
term in Eq. @) is not enough anymore. 

As discussed around Eq. |l3| ), we have the additional 
constraint that the map must preserve the value of To 
in order for the NLOi calculation to be consistent with 
the To resummation. This property was also important 
for the 1/2-jet separation to preserve the To resumma¬ 
tion. 

Technically, the only real-emission integral that we 
have to perform as part of the NLOi calculation is the 
integral in Eq. (Jl7). We perform this integral numeri¬ 
cally by Monte-Carlo integration over the radiation phase 
space for a given $1 event. The B 2 and C 2 terms in 
Eq. ( fl7| ) are allowed to use different maps, since all IR- 
safe maps agree in the IR-singular limit. However, we 
have to be able to invert both maps, i.e., we have to re¬ 
construct all $2 points that would project to the given 
$1 point. We then parametrize the radiation phase space 
for both maps in terms of a common set of radiation vari¬ 
ables, which are then Monte-Carlo sampled. While the 


The map itself is constructed as follows: We first clus¬ 
ter (by simply adding/subtracting four-momenta) the 
two partons in d> 2 that set the Ti value, i.e., that are 
closest according to the iV-jettiness metric. If the two 
final-state partons are clustered (FSR clustering), then 
this yields a $1 point with an off-shell final-state mo¬ 
mentum. If one of the final-state partons in $2 is closest 
to a beam direction (ISR clustering), then this yields a 
point with one of the incoming momenta off shell (and 
off axis). Next, we compute To on this off-shell $1 point. 
Finally, we construct an on-shell $1 point that has the 
same To value and the same of the vector boson. With 
g M and To given, the resulting point is fully determined 
by four-momentum conservation and on-shell conditions. 
Hence, by construction, the map also preserves the mo¬ 
mentum of the vector boson g M ($ 2 ) = ($ 2 ))- It is 

easy to see that this map is IR safe, and it is sufficiently 
simple that it can be inverted. 

For all ISR clusterings as well as FSR clusterings with 
both partons in the same hemisphere (defined by the 
vector-boson rapidity) this definition of 7))™ is identi¬ 
cal to the exact value of Toi^z)- The only case where 
it differs is for FSR clusterings where the two final-state 
partons are in opposite hemispheres. In the 7o-singular 
limit, this can only happen for soft emissions, and hence 
will in principle affect the constant terms in the two-loop 
soft function. (Collinear emissions in the beam functions 
always correspond to ISR clusterings.) However, this can 
only happen in a small region of phase space where the 
soft emitter was already accidentally close in rapidity to 
the vector boson. We therefore expect this to be a small 
effect, which we have verified numerically. We can there¬ 
fore safely ignore it for our resummation and simply take 
it into account as part of the nonsingular matching cor¬ 
rections. 

Finally, we need to comment on the ©^(d^) appearing 
in Eq. (12). It is easy to see that an on-shell massless $1 
point must satisfy To < |<fr|- However, there can be 
$2 points that can have a larger To value, and hence 
such <I> 2 points cannot be projected onto a massless <f>! 
point. The 0 7 "( < f> 2 ) is defined to be 1 for points that 
can be projected and 0 otherwise. Equivalently, these 
points also cannot be reached by the inverse map from 
splitting a *!>! point, and so are never part of the integral 
in Eq. ©■ Since the map is IR safe, the nonprojectable 
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$2 region is nonsingular and is included in Eq. (18) by 
the 1 — 0 7 ”($ 2 ) constraint. 


3. INTERFACING WITH A PARTON SHOWER 


which we define here in terms of the derivative of the 
Sudakov form factor and the normalized splitting proba¬ 
bility given in Sec. |2 B 4| 


From Eq. (42) and using the analog of Eq. (20) it is 


clear that the parton shower is unitary. Defining 


In the previous section we have given all required for¬ 
mulas for the jet cross sections dciQ 0 , dcrf c and dcr> 2 - 
As discussed in detail there, these jet cross sections in¬ 
clude the contributions of higher multiplicity phase space 
points, as long as the jet resolution variable for these 
points in phase space is below the 7jf ut value. In Table [i] 
we summarize how the phase-space points of different 
multiplicities contribute to the given jet cross section. 
Only the first three columns were necessary for the par- 
tonic calculation of the previous section. The last column 
shows instead the the constraints to be imposed on events 
with higher multiplicities. We will detail how these affect 
the development of the parton shower below. 

The purpose of the parton shower is to make the calcu¬ 
lation fully differential in the higher multiplicities, which 
can be viewed as filling the exclusive 0- and 1-jet bins 
with radiation, as well as adding extra jets to the inclu¬ 
sive 2-jet multiplicity. A parton shower acts on a par- 
tonic event by adding extra radiation in a recursive and 
unitary fashion. If the evolution variable of the parton 
shower were chosen to be iV-jettiness Q a single emission 
would be given by 


[dcq/v] = der/v 


d<f> 


tMT;T x , A) + £4(7^ ax , T n ) (42) 
) 0(7)r x > Tn > A)] , 


d'b./v+i , 

x ~rv- V($N+1 


N 


where T^ ax ~ 0{Tn-i) is the maximum value of Tn 
that can be reached given the value Tn-i that charac¬ 
terizes the hardness of the configuration before the emis¬ 
sion and A is the shower cutoff. This means that the 
parton shower either keeps the JV-body event unchanged 
with the no-emission probability determined by the Su¬ 
dakov factor Lbv(7)y aax , A), or it adds an extra emission 
according to the probability 

U' N (J^ ax ,T N )V^ N+1 ), (43) 


1 Here and in the following we treat Tn as the jet resolution scale 
of a single emission. Since Tn is a global variable that depends 
on all emissions, this warrants some explanation. What we in¬ 
tend here proceeds from the condition that for massless particles 
Tn(<&n+i) is always the total plus component of the closest pair 
of particles, using e.g. the metric defined in Eq. < |23| ), relative to 
the direction of the sum of their momenta. Moreover, 7 ~n —>■ 0 
when <&n +l —y Near the singular limit, where the shower is 

important and a good approximation of the underlying physics, 
one can thus assume that the direction of the axes entering the 
definition of Tn is aligned to the direction of the N hard par- 
tons. In this limit, a singular emission does not change these 
directions. Therefore, the value of Tn does in fact represent the 
hardness of the emission and can be used as evolution variable. 


da£ +1 = da n U' N (T^ x ,T N )nd> N+ 1) 

x 9(T^ ax > Tn > A), (44) 


d<f> 


N 


the recursive shower can now be written as 

S[da N ] = dexjv U N (V? ax , A) + S[daf, +1 \ . (45) 

Note that the definition of the parton shower requires 


a phase-space map. This is clear from Eq. (441, which 
depends on dd’Ar+i/d'hjv- 

The above equation assumes that the evolution vari¬ 
able of the parton shower is V-jettiness. Pythia 8, which 
we interface to for showering, as well as any other par- 
ton showers currently existing, use a different evolution 
variables, e.g. the transverse momentum or the angular 
separation of the emission. However, one can imagine 
taking the output of any of these showers and recluster 
the partons according to the iV-jettiness metric defined 
in Eq. p3| . The resulting history of splittings is equal 
at LL with the recursive expression given in Eq. (451. 
In the remaining discussion, we will work with Eqs. (42) 


and (45), and it is understood that this is valid only after 


the aforementioned reclustering. Note that this construc¬ 
tion implies a recursive definition of iV-jettiness, which in 
general is no t identical to the standard definition given 
It does agree, however, with the 7jf R defi- 


Sec. 


2B1 


nition discussed in Sec. 2 B 6 which is used in Geneva. 


The parton shower is not allowed to affect results of 
the jet cross sections at the accuracy they were calcu¬ 
lated in the previous section. Starting from this condi¬ 
tion, let us first consider the constraints we would have 
to impose on a parton shower that is strongly ordered 
in iV-jettiness, as discussed above. Strong ordering in 
Tn means that the Tn value at each emission is much 
smaller than the T^ ax upper limit determined by the 
Tn-i value of the previous emission. In this case, we 
would need to ensure that at each stage in the shower 
each partonic phase space point satisfies the constraints 
given in Table [T] For example, when showering a 0-jet 
event defined by <I>o, the phase space point reached 
after the first emission must have 7o(<E > i) < 7jj cut and 
must be projectable onto $ 0 - Each subsequent emission 
only needs to satisfy To($n) < 7jj cut . When showering 
a 1-jet event defined by $i, the $2 point obtained after 
the first emission would need to have 7i($ 2 ) < Tf ut , and 
$ 2 would have to be projectable onto with the 70- 
preserving map discussed in Sec. |2 B 6| This implies that 
the d >2 point obtained after one emission has to satisfy 
7o(4> 2 ) = 7o(4>i) as well as ^($ 2 ) = g M ($i), be¬ 
ing the vector boson four-momentum. Each subsequent 
emission then only needs to satisfy 7i($jv) < Tf ut . 
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4>o 

4>i 



do#7d$ 0 

All 

6»r 0 (4>i) and 0 map (4> o ; 4>i) 

6 , r 0 (4>2) 

9t 0 {$n) 

dcrfVdSq 

- 

0r o (4>i) or0 map ($i) 

9t 0 ($ 2 ) and 9 Tl ($ 2 ) and0 map (4>i; $ 2 ) 

9t 0 ($n) and 9r t ($jv) 

dcr| C 2 /d4«2 

- 

- 

e To ($ 2 )and [ 0 Tl ($ 2 )or 0 map ($ 2 )] 

9t 0 ($ n) and 9r x (4*jv) 


TABLE I. Table showing how the different partonic phase-space points are contributing to the different jet multiplic¬ 
ities. Here we have defined 9 Tn ($m) 3 9[Tn{$m) < 7^ ut ], 9 Tn {$m) = 9[Tn{$m) > U ut ], 0 map _(4«o; $1) = 
[4>i projects onto <f> 0 ], 0 map (4>i; 4> 2 ) = [$2 projects onto $ 1 ], 0 map ($i) = [<E>i does not project onto any $ 0 ], #m ap (4>2) = 
[$2 does not project onto any<f?i]. 


For a shower ordered in a different evolution variable, 
very similar restrictions hold, however the special con¬ 
ditions regarding the phase space maps are not required 
on the first emission, but rather on the emission with 
the largest value of the jet resolution scale, which can 
happen at a much later stage in the shower. Failing to 
properly account for this could affect the accuracy at 
which observables are predicted and ultimately spoil the 
leading-logarithmic terms and the color-coherence effects 
that are built into the shower Sudakov factors. 

A solution to this problem is to ensure that the first 
emission has indeed the largest value of the jet resolu¬ 
tion scale. This can be achieved by performing the first 
emission in Geneva, with a probability obtained from a 
simple LL Sudakov factor and using a splitting that pre¬ 
serves the map that was so carefully constructed. One 
would then only need to ensure that subsequent emis¬ 
sions from the parton shower have a lower jet resolution 
scale than the first, which can be achieved by properly 
setting the starting scale for the showering and applying 
standard veto techniques. 


The general idea of this approach is to split an iV-jet 
differential cross section into two pieces 


1 “C (]rr MC 

^-(TT;Aw) = ^:(n ut )U N (T^,A N ), (46) 


d4> 


N 


N 


da^N+i 

d$jv+i 


(t n > A N ,n ut ) 


d Tn 


<\rr M0 

aa N^-N 


d$ 


N 


(TT;Tw) 


P($N+l) 


x e(n ut >Tn> A N ) : 


(47) 


using a Sudakov factor Un and splitting probability V to 
define a branching and no-branching probability as in a 
unitary parton shower. The first term is now an exclusive 
iV-jet cross-section with a jet resolution at Ajv, which can 
be much smaller than the T™ 1 value used in Geneva. 
This process can also be iterated further. 

Applying this to our 4> 0 and 4>i events, and adding all 
emissions up to two partons, we find 


dag 10 

d'&o 


/' 7 “cut. 

V '0 5 


Aq) 


daf° 

d<f>i 


(T 0 >A 0 ,7? ut 


T 'CUt 
1 5 


Ai) 


drr MC 

ua >2 

d$2 


(To>A 0 ,7j>A 1 ,7? ut ,'7r t ) 


^(r o cut )t/o(T 0 cut ,A 0 ), 

^(T 0 > rrsrr*) C/i(71 cut , AO 

d d<T MC 

+ Tr 7ff|-(' 7 ? ut ; To) V(*i) 0(T o cut > % > A 0 ) t/i(77 iax , A,), 
a /0 

d (j MC 

-^(T 0 >r 0 cut ,ri >r 1 cut ) 

d d it m0 ~ 

+ cm ik {To > Ao ’ r ° cut ’ r i ut ’ T ^ v ^ 7T“ > Ti > Ao. 


(48) 


(49) 


(50) 


The 7j max in the theta-function in the last line above is 
either 7i cut or 7j max , depending on whether the derivative 
is on the first or the second term of Eq. (49). Once the 
first emissions have been performed as described above, 
the only remaining constraint on the parton shower is 
what is given in the last column of Table [TJ namely the 
constraint on the jet resolution variable. In practice, we 


will only apply this to the $1 events, since as discussed 
below, the showering of the $0 events is not affecting the 
accuracy of the 7o distribution. 

In the remainder of this section we give an argument 
why the parton shower acting on the jet events as de¬ 
scribed above does not affect the NNLL'+NNLO accu¬ 
racy. The $0 events include all phase-space points for 
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which To < To ut , and by definition Geneva is only pre¬ 
dicting the normalization, and not the distribution in To- 
Therefore, as long as we constrain the parton shower to 
only shower these events up to 7^ cut , the shower will ex¬ 
actly provide the missing To shape below 7^ cut which is 
integrated over by Geneva. On the other hand, due to 
the unitarity constraint, the shower will not change the 
normalization of the events below 7 J ut , which is thus still 
the NNLL'+NNLOo one predicted by Geneva. 

By performing the first emission of our $1 events as 
described above, and taking A-| as small as possible (we 
use a value of 0.1 GeV), the remaining 1-jet cross section 
is only about 0.1% of the total cross section. The theo¬ 
retical problem of how the shower affects the accuracy of 
this small contribution to the total cross section, which 
requires imposing all the constraints detailed above, can 
then be ignored for practical purposes. 

What is therefore left to show is that the shower does 
not affect the To distribution when showering $2 events. 
The To distribution of the original 2 -body phase space 
point is given by 

j| = /" 2 5S S|r “- 75(3,2)1 • (51) 

while the expression after the first emission done by the 
shower is given by 

^ = ^ M 6 l T ° ~ W 2 )] 

+ | d $ 3 ^^ro-r 0 ( $3 ) ]. ( 52 ) 

Clearly, thanks to the infrared safety of To, in the soft 
and collinear limit of $3 one has 


of this section. The reader can easily convince themselves 
that the argument given below will not be affected by this 
approximation. 

This allows us to write 

^ = / d<f>2 £ C/2 ( T “ aX ’ A2 ) «^) - 
r a s 

+ J d$ 3 ^«$2)-To + aT 2 ($3)]. (56) 

Thanks to the normalization of the splitting probabili¬ 
ties, one can easily perform the integrals over the two 
radiation variables other than T 2 to obtain 


do 5 

dTo 


= / d<f>2 


dcr 2 

(IT; 


C2(T 2 max ,A 2 )<5[To($2)-To] 


dcr 2 


r 1 


+ J d *‘^L d7 ^ (7?, “ r2) 

x«[7i(®j)-7; + a7y. (57) 

Next, the delta function is Taylor expanded to obtain 


da _ da s 

dTo~dT> 


= a d<E 


da 2 , 


S'IT 0 -T 0 (d> 2 )] 


d < f > 2 

rTl dTTU^TrTT) 


' a 2 


= a J dTo' ^ (T 2 ) (To') 5’[To - T 0 '] 


dT 0 


'0 
da 
d To 


(T 2 ) (To) 


(58) 


where we have defined the average T 2 value 


70(^3) = 7o( < i > 2) • 


(53) 


If the phase-space map of the parton shower were to pre¬ 
serve the value of To, the relation Eq. (53) would remain 
exact even away from the collinear and soft limits. The 
first emission done by the shower would then not change 
the To spectrum 


da s dcr 
d7o~dTo’ 


(54) 


If the phase-space map of the shower does not pre¬ 
serve To, instead, a difference between these two values 
exists. We can assume this difference to be proportional 
to T( < h 3 ), which measures the distance from the soft and 
collinear limits 


To($ 2 ) - To($ 3 ) =a($3)T 2 ($ 3 )- (55) 

Here a(<f> 3 ) is well behaved in the singular limit of 
$ 3 , because all the singular behavior is incorporated in 
7 2 ( 4 > 3 ) -A 0. For the sake of simplicity and to avoid 
rendering the notation too heavy we will drop the ($ 3 ) 
dependence of a and treat it as constant in the remainder 


C ! O. 


(T 2 ) = 


dT 2 T 2 U' 2 (TT^,T 2 ). 


(59) 


In the last lines of Eq. (58), we have written (T 2 ) (To) to 
make the dependence of the average value of 2 -jettiness 
on To explicit. In fact, one can derive this dependence 
using a simple LL Sudakov factor for the parton shower. 
Let us first consider a single emission governed by 


f/(T max ,T) ~ exp 


-C * In 2 - 
7r V 


T 


(60) 


with C the appropriate color factor. We can now easily 
show that 


xy—max 

(T) = f d TTU'(T max ,T) 

Jo 


T n 


e ir/(4a B C) 7r g r f c 


1 - 


\2^C J 


2 


C a s 

2 - -T 

7 r 


TO(al). 


(61) 
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Iterating this over two emissions, we find 

<r 2 )= [ T dTiU[(r 0 ,Ti ) Attu^t) 


C 2 a s 


/ViTiC/^TcTi) 

Jo 


C l( 7 2 a s 2 3 

' 4-^— T + CK^s) > 


(62) 


where we used the scaling 7)(p ax ~ Tn-i- 

We can also numerically extract the dependence of (75) 
on To from an actual shower in the following way: Start¬ 
ing from the showered expression d erf /d<f> 3 defined as in 
Eq. (44), we have 


/d$3^f<5(To($ 2 )-T 0 )75 

J d$3 


clcr 2 


r •2 


= /d$ 2 —^5[To-r 0 ($ 2 )] / dT a T a 05(77 

7 Ct® 2 7 a 2 

= JdT^6[To-Tj} (75) (T') 

= P o (mr o) . 


,75) 


(63) 


This gives 

(75) = 


dcr 2 

d7o 


-1 


« i 5 

d< i>3^^(ro($ 2 )-ro)r 2 . (64) 


Taking the $ 2 events from Geneva, running them 
through Pythia 8 , and calculating this ratio for vari¬ 
ous values of To, we obtain a very good fit to a straight 
line with 


(75) « 0.06 T 0 . 


(65) 


Taking advantage of the linear scaling of (75) with 
To, and dropping the explicit dependence in the notation 
from here forward, we can now rewrite Eq. (58) as 


der 

d7o 

This leads to 


dof 

dT 0 


= —a 


d da 

d% ° dT 0 


(75) 

To 


da _ d a s 

dTo dTo 

da 

dTo 


= /(To) 


(75) 

% 


■ ( 66 ) 


(67) 


where 


f(T 0 ) = -a- 


In 


dcr 

din To 


( 68 ) 


din To 

Since the dominant contribution to dcr/dln7o is given by 
dcr 


—a s ln7o e “ sln r ° 


(69) 


we find 


/(To) 


In To ' 


(70) 


This gives our final result for the change of the To spec¬ 
trum after the first emission done by the parton shower 


da 

dTo 


dcr 

dTo 


1 (75) a 2 


dr. ln7o T InT 


(71) 


Since the parton shower is strongly ordered in Tn, in 
the sense discussed at the beginning of this section, all 
the subsequent emissions will not alter this result to the 
order we are working. 

This result can be compared to the relative uncertainty 
in our NNLL' calculation of the To spectrum, which can 
be estimated considering the dominant term beyond the 
NNLL'order. This term scales as a^/To- Comparing this 
to the first dominant one that we include in the NNLL' 
prediction, we find the relative uncertainty to be 


A 


dcr 1 ' 


dTo 


a s/To 


da NNLL ' 

dTo 


c*s In To/To In To ' 


(72) 


Upon inspection of Eq. (|7T|) and Eq. (72), we can see that 


the shift in the spectrum due to Pythia’s showering is 
of the same size as the missing higher-order terms in the 
perturbative calculation. Therefore, the showering only 
affects our prediction beyond the claimed accuracy. 


4. COMPARISON WITH DEDICATED 
PERTURBATIVE CALCULATIONS 

In this section, we present the results of Geneva 
and the comparisons with dedicated FO and resummed 
predictions. All calculations are performed for pp col¬ 
lisions with 7TeV center-of-mass energy. We use the 
CTIOnnlo [73] set from Lhapdf6 J75] as our default 
PDFs, together with its default value of a s (Mz) = 0.118. 
We include the full interference effects between Z/T and 
restrict the invariant mass for the dilepton pair to the 
range 60 < < 120 GeV during the generation of 

events. The other parameters relevant for our calculation 
are 

M z = 91.1876 GeV, T z = 2.4952 GeV, (73) 
sin 2 Off = 0.2226459, a~^(M z ) = 132.338 . 

The Z couplings are given by 


e 

9e/q = sin 9^ cos 9$ 


Q) 


- Qe/q sin 2 6w 


(74) 


where i/q denotes the given left or right component of a 
lepton or a quark and e = \J 47r a em (M z ). For the par- 
ton shower, we use Pythia 8.205 [531155] . with a slightly 
modified version of Tune 4C, namely with a s = 0.118 


din To 
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FIG. 3. Validation of the To spectrum in Geneva. The analytic NNLL'+NNLO To resummation is compared to the partonic 
and showered results of Geneva, but before addition of nonperturbative effects. 



pp —> Zj'i —> e + e (7 TeV) 
Geneva+Pythia8 NNLL^+NNLO 0 

- Showered 

Hadronized 




To [GeV] 


FIG. 4. Comparison of the showered and hadronized To spectra in Geneva. 


for both the ISR and FSR shower, p = 4 GeV, and 
the hard primordial kr = 2.5 GeV. Unless otherwise 
stated, the Geneva predictions include showering and 
hadronization effects, but presently do not include mul¬ 
tiple parton interactions. 

We first present the To spectrum, and show that we ex¬ 
actly reproduce the input resummed spectrum, and that 
the parton shower does not modify this distribution as 
required. 

To validate the FO accuracy of Geneva, we next con¬ 
sider the rapidity distribution of the iGboson, as well 
as that of the negatively charged lepton from the Z de¬ 
cay. These distributions are known to NNLO accuracy. 
As discussed in Sec. [ 2 J Geneva reproduces all 0(a%) 
contributions up to small power-suppressed effects, and 
these observables allow us to test the size of these effects. 
In addition, it allows us to validate that the procedure 
described in Sec. 12 B 21 to turn off the resummation and 
numerically reproduce the total cross section from the 
integrated resummed spectrum works. 

Finally, to study the accuracy of our resummation, we 
show the results for various observables that are sensitive 
to the higher logarithms included in Geneva. The first 
is the 0-jet cross section, which only contains events with 
jets below a given pr cut and is known analytically to 
NNLL. The corresponding observable in Higgs produc¬ 
tion is an important ingredient in Higgs analyses. The 
other two resummation sensitive observables we consider 


are the transverse momentum distribution of the Z bo¬ 
son, as well as the <f>* distribution [76]. For all these 
observables we compare our results to the analytic re¬ 
summed results at NNLL matched to NNLOq or NLOi. 


A. The To Spectrum 

The 0-jet resolution variable chosen in Geneva is 0- 
jettiness, such that partonic configurations with 7o < 
7jj cut are part of the exclusive 0-jet MC cross section, 
while partonic configurations with 7o > 7o cut contribute 
to the inclusive 1-jet cross section. Since the dependence 
on the 0-jet resolution variable at NNLL' accuracy is used 
as input in Geneva, it must reproduce the To spectrum 
at NNLL' accuracy both for the central value and for the 
event-by-event scale uncertainties. 

In Fig. [3] we show the To spectrum in the peak, tran¬ 
sition, and tail regions. All the curves are cut off at the 
To ut = 1 GeV value used in the simulation. We show the 
analytic expressions with a red band, where the uncer¬ 
tainties are obtained using the same profile scale varia¬ 
tions discussed in Sec. |2 B 2| The green histogram with 
error bars shows the results of Geneva before attaching 
the parton shower. One can see that Geneva is repro¬ 
ducing the analytic results exactly, as it should, for both 
the central value and the scale variations. Finally, the 
blue histogram shows the results of Geneva after the 
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FIG. 5. Comparison of Geneva with the NNLO rapidity 
distribution of the vector boson. The orange curve shows the 
results from Dynnlo, while the black histogram shows the 
Geneva result. For Geneva, the uncertainties shown are the 
FO uncertainties as described in the text. 


parton shower. As required, the parton shower changes 
the % spectrum only within the perturbative uncertain¬ 
ties. In fact, over most of the To range the spectrum is 
essentially unchanged. 

The nonperturbative effects are shown in Fig. [4] where 
we compare the predictions for the To spectrum after the 
shower (again in blue) and after the addition of proton 
remnants, intrinsic kx smearing, and hadronization by 
Pythia 8 in orange. These effects change the To distri¬ 
bution significantly in the peak region, while they become 
power corrections in the transition and tail regions. This 
behaviour is precisely as dictated by factorization, from 
which one expects that these effects should behave for To 
analogous to thrust in e + e _ . Comparing to our e + e _ re¬ 
sults gn, this is precisely what we observe. A benefit of 
the Geneva approach is that it allows to directly com¬ 
bine the higher-order analytic resummation with these 
nonperturbative corrections provided by the hadroniza¬ 
tion model in Pythia 8. 


B. Partonic NNLOo Observables 

We now show that Geneva reproduces fully inclusive 
observables at NNLO accuracy, by comparing to dedi¬ 
cated NNLO calculations. In this section, we only con¬ 
sider the profile scale variations that reproduce the FO 



FIG. 6. Comparison of Geneva with the NNLO rapidity dis¬ 
tribution of the negatively charged lepton. The orange curve 
shows the results from Dynnlo, while the black histogram 
shows the Geneva result. For Geneva, the uncertainties 
shown are the FO uncertainties as described in the text. 


scale variations, as described in Sec. |2 B 2| 

In Fig. [5] and Fig. [6] we show the result for the rapid¬ 
ity distribution of the vector boson and the negatively 
charged lepton from its decay, respectively. The orange 
band shows the NNLO result from Dynnlo [77] . We 
show the results of Geneva as a black histogram, with 
the error bars representing FO uncertainties as described 
above. In the lower part of each plot, we show the ratio 
to the Dynnlo central value. 

The central value of Geneva agrees very well with 
the fixed NNLO prediction. The perturbative uncertain¬ 
ties from Geneva are also in reasonable agreement with 
those from Dynnlo. The few fluctuations observed in 
the plot are of statistical nature, as evidenced by the fact 
that they grow larger toward larger rapidities, where the 
statistics is poorer. The rapidity distribution of the vec¬ 
tor boson has also been validated against the independent 
NNLO calculation provided by Vrap m- 

In Fig. [7] we show the result of the transverse momen¬ 
tum distribution of the negatively charged lepton. For 
PTi < rrif+g-/ 2 this observable is a true NNLO distri¬ 
bution and Geneva agrees very well with the NNLO 
prediction. The region above and below m^+g -/2 is very 
sensitive to Sudakov shoulder logarithms E23- It is known 
that the FO calculations perform very poorly in this re¬ 
gion and fail to provide an accurate description of the 
physics. On the other hand, Geneva will have some of 
these logarithms resummed and should therefore provide 
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FIG. 7. Comparison of Geneva with NNLO transverse- 
momentum distribution of the negatively charged lepton. The 
orange curve shows the results from Dynnlo, while the black 
histogram shows the Geneva result. For Geneva, the uncer¬ 
tainties shown are the FO uncertainties as described in the 
text. 


a more accurate prediction. Finally, the region above 
Ptz > to^+£-/ 2 can only be populated if the transverse 
momentum of the vector boson is nonzero, making this 
region only NLO accurate. Close to pre ^ m ^+^-/2 
Geneva lies above the FO prediction and converges back 
to the FO result for large values of pn- This is likely still 
the effect from the Sudakov shoulder at m^c- / 2. Since 
the cross section above is much lower, any spillover from 
below TO|+£-/2 caused by the resummation can have a 
large relative effect. 


C. 0-Jet Observables with Resummation 

Since we resum our 0-jet resolution variable 7o to 
NNLL', it is interesting to study how accurately Geneva 
is able to predict other observables sensitive to the 
0/1—jet separation. Since the distribution for other 
0—jet resolution variables is only obtained indirectly, via 
the 7o spectrum made fully differential by using split¬ 
ting functions and a parton shower that are only correct 
to LL, we cannot formally claim the same NNLL' accu¬ 
racy for these observables. However, since the overall 
distribution of events in the 0/1-jet region of phase space 
is clearly improved, we expect some of this accuracy to 
carry over to other observables as well. Hence, we expect 
to get predictions that are numerically much closer to 



Qt [GeV] 


FIG. 8. Comparison of Geneva with the qr distribution from 
DYqT. The NNLL+NLOi (NLL+LOi) results from DYqT 
are shown in blue (green), and the Geneva results are shown 
in black. The Geneva results are much closer to DYqT’s 
NNLL results than NLL results. 

NNLL resummed results also for other 0-jet observables. 
The same behaviour was already observed for e + e _ 24] . 

In Fig. [8] we show the transverse momentum distribu¬ 
tion of the vector boson compared to its analytic NNLL 
resummation from DYqT. The predictions of DYqT 
have been manually switched to agree with the FO re¬ 
sults in the tail, according to the procedure outlined in 
Ref. [S3- We see that Geneva agrees reasonably well 
within the perturbative uncertainties with the NNLL re¬ 
summed resultj^] The NLL result from DYqT has a sig¬ 
nificantly different shape, and the Geneva prediction is 
certainly in much better agreement with the NNLL pre¬ 
dictions than the NLL ones. 

A similar prediction for the transverse momentum dis¬ 
tribution, but fully differential on the momenta of vector 
boson decay products has been presented in Ref. m • 
This allows the direct comparison including the accep- 


2 Although Geneva’s perturbative uncertainties appear smaller 
than DYqT’s at very low qT , this should not be misinterpreted as 
Geneva being more accurate here. This is mostly an artifact of 
smearing out the uncertainties from a range of 7o values, which 
tends to reduce their numerical impact. In addition, we have not 
included here any uncertainties associated with the interface to 
the parton shower and the showering itself. Such uncertainties 
should eventually also be included as an additional part of the 
(resummation) uncertainties. As this will require a dedicated 
study, we leave this for future work. 
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FIG. 9. Comparison of Geneva with the analytic qr distribu¬ 
tion at NNLL+NLOi from Ref. [50]. The analytic results are 
shown in blue, and the Geneva results are shown in black. 



FIG. 10. Comparison of Geneva with the analytic cp* distri¬ 
bution at NNLL+NLOi from Ref. ;80]. The analytic results 
are shown in blue, and the Geneva results are shown in black. 


tance cuts used by experimental analyses. Using the 
same lepton cuts as in that study, we show our compari- 


FIG. 11. Comparison of Geneva with the 0-jet cross section 
as a function of p from JetVHeto [HT| at NNLL+NNLO. 
The analytic results are shown in blue, and the Geneva re¬ 
sults are shown in black. For Geneva, the uncertainties are 
the FO uncertainties only, see text for details. 


son in Fig. [9] Again, we observe a fairly good agreement 
with the analytic NNLL prediction matched to NLOi. 
Another variable, quite similar to the transverse momen¬ 
tum of the vector boson, is the cp* between the two lep¬ 
tons, with the precise definition of cp* given in [76]. The 
comparison of Geneva to the NNLL+NLOi calculation 


of Ref. [80] is shown in Fig. 10 and we again observe 
good agreement. 

Finally, we show the result for the exclusive 0-jet 
cross section as a function of p™ 1 in Fig. 11 where 
the 0-jet sample is defined as all events containing no 
jets with pt > pjP. The jets are reconstructed with 
the anti— kp algorithm [33 J as implemented in Fast- 
JET [5U[55] , within a radius R = 0.4. We find good agree¬ 
ment between Geneva and the dedicated NNLL+NNLO 
calculation given by JetVHeto [ST] within the pertur¬ 
bative uncertainties. For this plot, we use the FO scale 
uncertainties discussed in Sec. |2 B 2[ such that the uncer¬ 
tainties at large p™ 1 are estimated correctly and thus pre¬ 
cisely reproduce those of JetVHeto. At small p™‘ they 
are now underestimated and here the resummation un¬ 
certainties should be added. The better agreement with 
the NNLL+NNLO prediction compared to the lower or¬ 
der NLL+NLO one, especially in the large p™ 1 region, 
is of course driven by the correct inclusion of the NNLO 
corrections in Geneva. 
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FIG. 12. Comparison of Geneva with the rapidity of the 
vector boson from the ATLAS study in [SB]. The Geneva 
results with default values of a s (Mz) and with a s (Mz) = 
0.1135 are shown in orange and green, respectively, while the 
ATLAS points are shown in black. 


5. COMPARISON TO LHC DATA 

In this section, we compare the results from Geneva 
with measurements using data collected during the 7 TeV 
LHC Run I. We use Geneva +Pythia 8 to produce fully 
showered and hadronized events in HepMC format, and 
then use the Rivet routines provided by the experimen¬ 
tal collaborations to produce the histograms that can be 
compared with the experimental measurements. By us¬ 
ing Rivet S3, we ensure that event selection cuts and 
object definitions agree with what was used in the ex¬ 
perimental analyses. We use the same run that was used 
for the comparison to the perturbative calculations, such 
that all parameters are identical to those given in Sec. [4] 

In addition, we also include in our data comparisons 
Geneva results with a lower value of a s (Mz ) = 0.1135, 
since as we will see, the value of a s (Mz) can have a no¬ 
ticeable effect on resummation-sensitive observables. In 
the following, the results with a s {Mz) = 0.118 will be 
shown by the orange bands, while those with a s (Mz) = 
0.1135 are shown by a green band. For the latter re¬ 
sults, we only change the a s (M z ) value, holding all other 
parameters including the PDF set fixed]^] This value 


3 We hasten to add that using a different a a {Mz ) than the de¬ 


is motivated by the fact, that all a s (Mz) determina¬ 
tions from LEP data that include higher-order resum¬ 
mation together with some form of analytic treatment of 
nonperturbative hadronization effects, find considerably 
lower values of a s (Mz) than 0.118 [57HT01 [58H5P] , see 
also Ref. [ST] for determinations from DIS. Consistent 
with these findings, in Ref. [21] it was observed that the 
higher-order resummation using a s (Mz) = 0.1135 cou¬ 
pled with Pythia 8’s hadronization model gave excellent 
agreement with data for e + e _ —> 2-jet resummation ob¬ 
servables, namely thrust, C-parameter, heavy-jet mass, 
and jet broadening. On the other hand, the predictions 
with a s (Mz) = 0.118 gave consistently spectra that were 
somewhat too hard, i.e., shifted from lower to higher val¬ 
ues compared to data. We will observe precisely the same 
qualitative trend here. 

The rapidity distribution of the vector boson was mea¬ 
sured by ATLAS jS6], CMS (92], and LHCb [93]. Since 
we already saw that Geneva agrees very well with the 
NNLO predictions, these comparisons are mostly driven 
by the PDFs we use. We show the comparison to the ra¬ 
pidity spectrum from ATLAS and CMS in Figs. [12] and 
[13] The normalized results from CMS are reproduced 
very well for the entire Y range. The ATLAS results are 
not normalized, and while Geneva agrees very well at 
low Y values, it slightly overshoots the data (by about 
2 standard deviations) for larger Y. A similar discrep¬ 
ancy was observed in Ref. [SB] for the MSTW08, HERA- 
PDF1.5, and ABKM09 PDF sets. The comparison with 
LHCb is shown in Fig. [14} The agreement at moder¬ 
ate values of Y is good here too. For forward events 
with Y > 3.25 our predictions show a trend to be lower 
than the data. This discrepancy was already noted in 
the LHCb paper [93] when comparing with the NNLO 
calculation provided by Fewz [91]. As we would expect, 
the lower a s (Mz) value makes almost no difference for 
the inclusive predictions. 

A measurement of the production cross section of a Z 
boson in association with jets was presented by ATLAS 
in Ref. [95] , for both inclusive and exclusive jet cross sec¬ 
tions. We compare our results with the measurements 
having up to 2 jets in Fig. [15] We choose to limit our¬ 
selves to up to 2 jets because any additional jet would 
only be provided at LL accuracy by Pythia 8. Our pre¬ 
dictions for the 0-jet cross section agrees well with data, 
not only for the inclusive cross section, but also for the 


fault value provided by the PDF set is not problematic as it may 
seem at first. In the resummation-dominated region, this only 
leads to a small inconsistency inside the beam functions, where 
the terms that are supposed to cancel the PDF’s dependence 
will now use a different o; s value leading to a small noncancel¬ 
lation in the coefficients of some single-logarithmic terms. It is 
safe to assume that this will not be numerically relevant com¬ 
pared to the dominant a s dependence coming from the double- 
logarithmic beam function terms, as well as the hard and soft 
function, and all the resummation kernels. All of these pieces 
are completely unrelated to the PDFs, and are in fact the same 
as for the corresponding T 2 resummation in e+e — collisions. 
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FIG. 13. Comparison of Geneva with the rapidity of the 
vector boson from the CMS study in [92] . The Geneva results 
with default values of a a (Mz) and with a a (Mz) = 0.1135 are 
shown in orange and green, respectively, while the CMS points 
are shown in black. 



FIG. 14. Comparison of Geneva with the rapidity of the 
vector boson from the LHCb study in 93]. The Geneva 
results with default values of a a (Mz) and with a a (Mz) = 
0.1135 are shown in orange and green, respectively, while the 
LHCb points are shown in black. 


exclusive one, where resummation plays a role. The un¬ 
certainties for our predictions of the inclusive Z + 1-jet 
cross section are larger, since they are only predicted to 
NLOi, though they can still benefit from the resumma¬ 
tion of the 0/1 jet boundary. The separation into ex¬ 
clusive Z + 1-jet and inclusive Z + 2-jet is only at LL 
accuracy, and the FO accuracy of the inclusive Z + 2- 
jet cross section is only LO 2 , with correspondingly larger 
perturbative uncertainties. Geneva agrees well with the 
data, somewhat better for lower a s (Mz). 

In Ref. [95] , ATLAS also presented a measurement of 
the transverse momentum distribution of the hardest jet, 


and the comparison to these results is shown in Fig. 16 


The predictions from Geneva are in good agreement 
with the data. Below p ] ^ x < mz the agreement is no¬ 
ticeably better for lower a s (Mz). Above pl^ > 300 GeV 
the predictions tend to be higher than the data, but still 
consistent within the larger uncertainties. This could be 
due to the fact that we use a renormalization scale of 
p = mg+e- , while at such large transverse momenta a 
better choice might be p = p ] ^ x . 

Next, we compare our predictions to the </>* between 
the leptons, and the transverse momentum of the vec¬ 
tor boson. Both of these observables are 0-jet resumma¬ 
tion variables, since at low values they are dominated by 
events without any hard emissions. Both of these distri¬ 
bution are quite sensitive to the choice of parameters of 


the parton shower and the nonperturbative model used 
in Pythia 8. Without higher-order perturbative cor¬ 
rections included, the MC tune will partially adjust the 
available parameters to mimic missing higher-order per¬ 
turbative effects. This implies that when using Pythia 8 
in conjunction with Geneva, which includes much more 
higher-order perturbative information, a retuning of the 
parameters becomes necessary. We stress that no at¬ 
tempt at a systematic retuning of Pythia 8 has been 
done for this work, and we expect that a dedicated tune 
of Pythia 8 together with Geneva, will improve the 
data agreement. 

The (f >* distribution has been measured by ATLAS in 
Ref. [SB] and LHCb in Ref. [93] ■ The measurement from 
ATLAS is a normalized spectrum, while LHCb quotes 
the unnormalized distribution. The comparisons with 
these two measurements are shown in Figs. [T7] and |18| 
Geneva agrees well with the results from LHCb, though 
the measurement has relatively large uncertainties. Com¬ 
paring with the much more precise results from ATLAS, 
Geneva predicts a wider distribution, such that our 
predictions are below the data in the peak region for 
<f>* < 0.1, and above the data for larger values of 4>*. 
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FIG. 15. Comparison of Geneva with data on inclusive and 
exclusive jet cross sections from the ATLAS study in [951 . 
The Geneva results with default values of a s (Mz) and with 
a s (Mz) = 0.1135 are shown in orange and green, respectively, 
while the ATLAS points are shown in black. 


FIG. 17. Comparison of Geneva with the (j >* distribution 
from the ATLAS study in [96]. The Geneva results with de¬ 
fault values of a s (M z ) and with a s (Mz) = 0.1135 are shown 
in orange and green, respectively, while the ATLAS points are 
shown in black. 



FIG. 16. Comparison of Geneva with the transverse mo¬ 
mentum of the leading jet from the ATLAS study in [951 . 
The Geneva results with default values of a s {Mz) and with 
a s (Mz) = 0.1135 are shown in orange and green, respectively, 
while the ATLAS points are shown in black. 


FIG. 18. Comparison of Geneva with the <j>* distribution 
from the LHCb study in m- The Geneva results with de¬ 
fault values of a s (Mz) and with a s (Mz) = 0.1135 are shown 
in orange and green, respectively, while the LHCb points are 
shown in black. 
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FIG. 19. Comparison of Geneva with the transverse mo¬ 
mentum of the vector boson from the CMS study in [92j . 
The Geneva results with default values of a s {Mz) and with 
a s (Mz) = 0.1135 are shown in orange and green, respectively, 
while the CMS points are shown in black. 


Discrepancies with this same trend and at a similar 
level were observed for other MC predictions in Ref. [5E]. 
The agreement with data is considerable better when 
choosing a lower value of a s (M z ) = 0.1135. 

As a final comparison we consider the transverse mo¬ 
mentum (qx) distribution of the Z boson. The qx dis¬ 
tribution was measured by CMS in Ref. [92] and by AT¬ 
LAS in Ref. [57] , where the latter also measured the qx 
spectrum in different bins of Y. The comparison of our 
results with these measurements is shown in Figs. 19 and 
[20| The situation is similar to the <j>* distribution dis¬ 
cussed above. For a s (M z ) = 0.118, Geneva is below the 
data in the peak region qx < 10 GeV, while it is above 
the data in the transition and FO regions for qx > 10 
GeV. This is of course not unexpected, since the observ¬ 
ables qx and cf>* are highly correlated. The same level 
of disagreement was already observed in Ref. m, where 
the data was compared with the NNLL+NLOi results of 
Ref. j8T|. As for the <jf distribution, the agreement be¬ 
tween Geneva and the data is noticeably improved for 
a lower a s (M z ) = 0.1135. 


6. CONCLUSIONS 


We have presented a combination of the fully- 
differential NNLO calculation for Drell-Yan production 


pp —>■ 7/Z —I £ + £~ combined with the NNLL' resumma¬ 
tion of 0-jettiness and interfaced with the parton shower 
provided by Pythia 8 within the Geneva Monte-Carlo 
framework. The starting point of Geneva is the formu¬ 
lation of Monte-Carlo cross sections, which are IR safe 
partonic jet cross sections, defined using a jet resolution 
variable 7jv, which is chosen to be V-jettiness in our im¬ 
plementation. We include 0-jet, 1-jet, and 2-jet cross sec¬ 
tions at NNLOo, NLOi and LO 2 . Furthermore, the de¬ 
pendence on the 0-jet resolution variable 7o is resummed 
to NNLL' accuracy, while that of the 1-jet resolution vari¬ 
able 7i is presently resummed to LL. 

Interfacing these partonic results with a parton shower 
such as Pythia 8 requires a careful treatment. First, one 
needs to deal with the fact that the parton shower evolu¬ 
tion variable is different from V-jettiness. Second, care 
has to be taken that the parton shower’s lower resum¬ 
mation accuracy and different phase-space map does not 
destroy the higher logarithmic accuracy of the calculated 
Monte-Carlo cross sections. We discussed in detail how 
to solve these issues by performing the first two shower 
emissions by hand using the jet-resolution parameter as 
ordering variable. We then showed that any subsequent 
showering by Pythia 8 does not affect the formal per¬ 
turbative accuracy included in Geneva. 

We have validated the FO perturbative accuracy by 
comparing our results for inclusive 0-jet observables to 
Dynnlo. We agree with this dedicated NNLO calcula¬ 
tion within the small perturbative (and statistical) uncer¬ 
tainties. We also studied how Geneva’s improved per¬ 
turbative accuracy in the resummation region from re¬ 
summing 0-jettiness to NNLL', translates to other 0-jet 
resummation variables, such as the transverse momen¬ 
tum of the vector boson qx , the </>* between the two lep¬ 
tons, and the exclusive 0-jet cross section as a function of 
the jet Since Geneva partially relies on the parton 
shower for these observables, they formally do not have 
NNLL' accuracy. Nevertheless, we find that Geneva re¬ 
produces dedicated NNLL resummations for these three 
observables rather well, and in particular the Geneva 
predictions are much closer to the exact NNLL+NLOi 
results than to the NLL+LOi results. This is a clear in¬ 
dication that the gain in resummation accuracy for 7o, 
when implemented into a fully exclusive prediction, also 
translates into more accurate predictions for other 0-jet 
resolution variables. 

Finally, we presented a comparison of Geneva with 
measurements by ATLAS, CMS, and LHCb using the 
7TeV LHC data. Since Geneva agrees with the NNLO 
results for the rapidity distribution of the vector boson, 
the agreement between Geneva and the LHC data is 
very similar to that observed when comparing to other 
NNLO calculations. Here, we note that experimental 
measurements of the lepton pxe spectrum would be very 
valuable. 

We find good agreement between Geneva and the AT¬ 
LAS measurement of exclusive jet cross sections and of 
the transverse momentum distribution of the hardest jet. 
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FIG. 20. Comparison of Geneva with the transverse momentum of the vector boson both inclusive (top-left) and in different 
bins of vector boson rapidity Y from the ATLAS study in [97]. The Geneva results with default values of a s (Mz) and with 
a s (Mz) = 0.1135 are shown in orange and green, respectively, while the ATLAS points are shown in black. 


For the transverse momentum of the vector boson or the 
(j>* between the leptons our predictions when using the 
default value of a s (Mz) = 0.118 are lower in the peak re¬ 
gion and higher in the transition region behind the peak 
compared to the measurements. The same trend has 
been observed before by other MC predictions and also 
dedicated resummed calculations for these observables. 


When using a lower value a s (Mz) = 0.1135, we observe 
a noticeable improvement in agreement with data for es¬ 
sentially all resummation-sensitive observables. A simi¬ 
lar effect was observed previously for a variety of e + e - 
2-jet resummation observables. We believe this deserves 
further attention. 

We encourage the experiments to also measure 7o-like 
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jet-based variables (see Ref. EDI) and to also perform 
differential measurements of resummation observables at 
much higher m^+g- values. Since the resummation type 
and regime for these cases are very different, they would 
add valuable complementary information for theory com¬ 
parisons. In general, continued precise measurements of 
resummation observables at the LHC are essential for our 
understanding of higher-order QCD effects in exclusive 
and differential observables, which ultimately will also 
provide important inputs for the theoretical interpreta¬ 
tion of Higgs measurements and new-physics searches. 
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